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l/-) ■ tively on the Dirichlet and Neumann boundary value problems. The equivalence of 

these methods is shown due to nontrivial duality relations between special boundary 
operators of these two problems. Previously known braneworld action algorithms in 
two-brane Randall-Sundrum model are generalized to curved branes with deSitter 
, and Anti-deSitter geometries. 
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1 Introduction 

Rapidly developing theory of braneworld phenomena [|1], 0, |3| requires new efficient 
methods of their description. These methods incorporate together with the well-known 
old formalisms, like the effective action approach, special new features associated with 
bulk/brane (boundary) ingredients characteristic of the braneworld scenarios. In this 
paper we want to focus our attention at the peculiarities of the Dirichlet and Neumann 
boundary conditions in braneworld setup, the way they arise in the course of calculating 
braneworld effective action and, in particular, at a sort of duality relation between these 
problems. This relation manifests itself in the equality of special boundary operators 
originating from Dirichlet and Neumann boundary value problems and serving as a kernel 
of the braneworld effective action in the approximation quadratic in fields. 

Braneworld effective action implicitly incorporates the dynamics of the fields in the 
bulk and explicitly features the boundary fields (being the functional of those) which, 
in braneworld scheme, are directly observable by the observer living on the brane. Such 
action, on the one hand, arises as a result of integrating out the bulk fields subject to 
boundary (brane) fields and, on the other hand, generates effective equations of motion for 
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the latter. This situation obviously suggests two strategies of calculating the braneworld 
action. One strategy consists in its direct calculation and, in the tree-level approximation, 
reduces to solving the equations of motion in the bulk subject to given boundary values 
- brane fields - and substituting the result into the fundamental bulk action. Thus, 
this strategy involves the Dirichlet problem. Another strategy consists in recovering the 
braneworld effective action from effective equations of motion for brane fields. Since the 
latter incorporate well known Israel junction conditions on branes [f|, this method relies 
on the Neumann boundary value problem. The consistency of these two methods, which 
obviously should lead to one and the same result, is far from being explicit. Here we show 
that they really match in view of a nontrivial relation between nonlocal brane operators 
arising as restriction of (properly differentiated) kernels of Dirichlet and Neumann Green's 
functions to the boundaries/branes. As a byproduct of the Dirichlet setup we derive the 
answer for the braneworld effective action in models with curved branes of deSitter and 
Anti-deSitter geometries. 

Curved brane models have been considered in many papers their very incomplete list 
consisting of []5|, [7|, P, PH. In contrast to the results of these works we concentrate 
on two-brane models with independent metric configurations on their branes. This is 



done along the lines of the previous papers [13], |14[ where the weak field (quadratic in 



perturbations) approximation for braneworld action was obtained on the background of 
flat branes. In distinction from 0, |l(| we obtain the nonlocal action nontrivially inter- 
twining metric fields on both branes. The interest in such a construction follows from the 
attempts of describing off-shell phenomena in braneworld physics, solving the hierarchy 
and cosmological constant problems, generating inflation in braneworld scenarios [[TJ |TB| ], 
etc. 

The paper is organized as follows. In Sect. 2 we define the braneworld effective action 
and sketch its calculation within the perturbation expansion with Dirichlet boundary 
conditions on branes. In Sect. 3 we apply this technique to the two-brane Randall- Sundrum 



model |12j to find its action in the approximation quadratic in the perturbations of metric 
fields on branes. This result generalizes the algorithms of the works |13|, [14] to curved 
branes of deSitter and Anti-deSitter geometries. Its derivation is based on the technique of 
linearized invariants of the diffeomorphism gauge transformations intensively used in the 
theory of cosmological perturbations [TR [L6|, ITTA. Sect. 4 gives a review of the alternative 



calculation |14] of the same result, based on its recovery from the effective equations 
of motion. In Sect. 5 we show the equality of the two manifestly different quadratic 
forms of the action by proving that their kernels are equal to one another - the main 
duality relation between specific operators of the Dirichlet and Neumann boundary value 
problems. Finally, in concluding section we discuss possible implications of the obtained 
results in braneworld scenarios with curved (inflating) branes. 

2 Braneworld effective action 

We will calculate the braneworld effective action defined by the procedure that was 



discussed in much detail in [14]]. In contrast to the conventional Kaluza-Klein reduction, 



usually used for the construction of the effective action [|1^, [Tj| , our action is a functional of 
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the fields that have direct physical and geometrical interpretation. They coincide with the 
restriction of the fundamental bulk fields $ to the branes, they are directly observable 
by observers living on these branes, and serve as boundary conditions for $. 

To be more precise, let $ be a set of fundamental dynamical fields on (d + 1)- 
dimensional spacetime Ai (bulk) with the full boundary dAA. In dAi we shall include 
not only the "outer" boundary (like asymptotic domains of the bulk spacetime) but also 
the branes - timelike <i-dimensional surfaces embedded into Ai and enumerated by the 
index i, dAi = [j^Aii. Boundary values of $ on dAii we denote by fa, $\ 9M . = 4>i, and 
call them brane (or induced) fields. In addition to the bulk (d+ l)-fields there are usually 
(i-dimensional fields tp which, by definition, live only on branes and which will be in what 
follows denoted as matter fields^. 

Thus, the full (d + l)-dimensional action is 

&[$,<p] = S (d+1) [$] + S mat [0,y?], (1) 

where the first term denotes the bulk part depending only on $, while the second term 
represents the matter action which depends only on induced (brane) <fi and matter fields 
ip. For brevity we omit the index i of brane and matter fields, so that (<ft, if) denotes the 
collection of fields (0j, <pi) associated with all the branes. 

Braneworld effective action is the result of integrating out the bulk fields subject to 
given values of brane fields 



exp (iS eff [0,y?]) = J D$ exp (i8[$,(p 



(2) 



Since the action of matter fields enters ([[]) additively and its arguments are not integrated 
over, it continues entering S [<f>, <p\- The rest of the latter is highly nontrivial, because it 
accumulates the result of the functional integration. When this integration is done within 
^-expansion the result reads 

S eff [0,y,] = S d \<P] + S loop [0] + S mat [0,y,], (3) 

where the tree-level part 

^[0] = S( d+1 )[$[0]] (4) 

is a result of substituting in the classical bulk action the solution $[0] of the following 
Dirichlet boundary value problem 

SS^ _ 

" °' (5) 

Km = <t>- 

Also, S"' 0015 ^] represents the loop part, of order h, which we shall disregard in what follows. 

1 Usually the role of bulk fields is played by the bulk spacetime metric, because string theory motivated 
braneworld framework does not admit the propagation of gauge (matter) fields in the bulk M . 
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When considering the full nonlinear theory ([!]) it is generically impossible to obtain 
the tree-level term (|4]) as an exact explicit functional of induced fields 0. Instead, one can 
develop a perturbation theory for (HI) in powers of the perturbation of boundary conditions 



+ h 



(6) 



on the background of some configuration 0° at which the effective action is exactly calcu- 
lable (or represents an irrelevant constant). This perturbation induces the perturbation 
H[h] of the solution of the boundary problem in the bulk (||) 



$[0] = $° + H[h' 



(7) 



where $° = $[0°] is the bulk background solving the classical equations of motion with 
background boundary conditions on branes = 0°, which is also supposed to be explicitly 
known. 

From (H) and (|7|) it follows that up to 0(/i 2 )-terms, the bulk perturbation H[(j)°,h] 
satisfies the following linear Dirichlet problem: 

FH = 0, 



H\ 



h. 



where F is the operator of small field disturbances in the theory with the action (Jl]). Since 
we assume that its Lagrangian contains at most first-order derivatives, F is a second-order 
differential operator with the kernel in the space of (d + l)-dimensional bulk coordinates 
X calculated on the background $° 



F5(X, X') 



(9) 



As any other second-order differential operator, F satisfies the Wronskian relation 
following from integration by parts in the equation 



(10) 



dV $i(F$ 2 ) - (F$i)$ : 



dS $i(iy$ 2 ) - 



M 



dM 



for arbitrary $ and for arbitrary test functions $i and $2- Here dV is the bulk 
integration measure and dS is the measure on the boundary, and we shall call W the 
Wronskian operator. It is of the first order in derivatives and, for regular situations, it 
necessarily contains the derivative normal to the boundary. Obviously, this operator is 
defined by the relation (|T0|) only up to addition of an arbitrary selfadjoint operator on 
the boundary. It can be specified uniquely by demanding that 



#*iF$ 2 = / dV®i{F$ 2 ) 



dS$i(W$ 2 ), 



(11) 



M 



M 



dM 



where the notation F implies that the derivatives of F are acting (in the sense of in- 
tegration by parts) so that the left-hand side represents the form bilinear in first order 
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derivatives. For $2 = $1 the integrand of the left-hand side is just the Lagrangian 
quadratic in the field and its first-order derivatives. In this case the operator W can also 
be obtained in terms of the canonical momentum conjugated to the field relative to "time" 
normally flowing to the boundary |HJ . 



This Wronskian operator allows one to formulate in a closed form the solution of the 
Dirichlet problem (j|) given the Green's function of this problem satisfying 

FG D (X,X')=S{X,X% 

(12) 

G D (X,X')\ xegM = 0. 
This solution reads as the following integral over the boundary dAi 

$(X) = - I dS'\G D (X,X')W'] h(x'). (13) 

J L \x'=x{x') 

dM 

Here X = X(x) denotes the embedding of the boundary dM. parameterized by internal 
coordinates x into the bulk with coordinates X, dS' = dS(x') is the surface integration 
element at x' and the arrow over W indicates that the derivatives of the primed Wronskian 
operator act to the left on the second (primed) argument of the Green's function. 
The expansion of the action in the field perturbation (0) begins with 



S d [(f>] 



S (d+i)[$o] + 1 J dVH*FH + 0[h s ], (14) 

M 



because the linear term identically vanishes in virtue of classical equations of motion for 
$° and, what is important, in virtue of the Israel matching conditions on branes, which 
annihilate the linear surface terms on boundaries even despite the fact that H\qm — h 7^ 0. 

Note that the quadratic part of this expansion is built with the aid of the operator F 
symmetrically acting on both fields - the consequence of the fact that the initial action 
did not contain second-order derivatives. Therefore, using the relation ( pi] ) and taking 
into account linear equations of motion for perturbations (§) one finds that the bulk 
part of the quadratic term vanishes and the quadratic approximation for the braneworld 
effective action reduces to the surface termQ 

S40] = S( d+1 )[$°]+i J dSH{WH). (15) 

dM 

Further, a substitution of (13) in this expression gives the general answer for the quadratic 
action 

S^[h] = -^ J dS J dS'h(x)WG D W(x,x')h(x'), (16) 

dM dM 



2 In what follows we restrict ourselves with the quadratic approximation and omit terms 0[h 3 } 
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where the kernel of the integral operation on the boundary WGj)W(a;,i') is given by an 
obvious sequence of operations with the kernel of the Dirichlet Green's function - acting 
on it from both sides with Wronskian operators and restricting the result to the boundary 
with respect to both arguments 



WG D W(x,x) = WG D (X,X')W 



X=X(x),X'=X(x>) 



(17) 



As mentioned above, for the case of several branes or spacetime boundary consisting 
of several pieces dM. = [j { dMi the boundary (brane) fields denote the collection (ft = (f> 1 . 
Correspondingly the brane embedding functions and respective perturbations acquire the 
same index, X = Xi(x), h = h l = h\sMi, an d the equations above imply the following 
obvious generalization 



dS — > J dS, 

dM % dMi 

WG c W(i, x) -> [WGflWj^i, x') = WG D (X, X')W 



(18) 
• (19) 



Under these replacements the expression 
matrix- valued kernel (|19|). 



X=X i (x),X'=X j (x') 

becomes a quadratic form with the nonlocal 



3 Two-brane Randall- Sundrum model 

Here we specify the construction of the above type on the example of the two-brane 



Randall- Sundrum model 121 with the action 



S[G>]= J d d+1 X VG (R(G) - 2A) - 2 J d d x^[K]-J2 J d d x^a t + 

BxZ 2 1 hi 1 b t 

+ J2f d d x^L™\g,^dcp). (20) 

The orbifold symmetry implies that the (d + l)-dimensional integration runs over two 
identical copies of the bulk B which is bounded by two <i-dimensional branes bj, % = ±, 
that can be regarded as the boundaries of B. Here G = Gab(X) (A = 0, 1, d) is the 
bulk metric and g = <?ig(x) (a = 0, 1, ...,d— 1) denotes the collection of induced metrics 
on % = ± branes, [K] is a jump of the trace of the extrinsic curvature on the brane defined 
as K = G V A n B, where is a covariant (d +l)-dimensional derivative and n is an 
inward pointing normal. With this definition the normals on two sides of the brane are 
oppositely oriented and the extrinsic curvature jump [K a p] actually equals the sum of 
the so-defined curvatures on both sides of the brane. The orbifold symmetry obviously 
implies J BxZ (...) = 2 J B (...) and [K] = 2K. Finally, A is the bulk cosmological constant 
and (jj are brane tensions. 

Our goal is to construct the perturbative braneworld effective action for such a system. 
The fundamental bulk field is the bulk metric $ = Gab(X), the role of brane field is played 
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by = g^g(x) - the induced metric on branes. The first three terms of ( p0|) comprise the 
bulk part of the action and the last term is the action of matter fields located on branes|. 

We choose as the background $° the well-known Randall-Sundrum solution in this 
model. In the setting generalized to the case of curved branes this solution arises as 
follows. Assume that there exists a coordinate system on the (d + l)-dimensional bulk 
X A = (x a , y) in which the background metric solving the vacuum Einstein equations has 
the form 

g%(x) = ( a ^y)^) (21) 

Two branes bounding the bulk manifold B are the hypersurfaces of constant y located 
at y = y + and y = so that B has a topology of bxl with x a - the coordinates on 
ci-dimensional b and I - the interval of the coordinate y G [y+, y_]. In the wording of the 
previous section the boundary conditions 0° for this bulk background $[0°] are two in- 
duced metrics conformally equivalent to one another g a „ = a 2 (y±)g a p(x) and conformally 
related to some fixed metric g a p{x). With the ansatz ([H]) this metric is restricted by the 
condition of constant scalar curvature. Indeed, from the ^-component of the vacuum 
Einstein equations in the bulk it follows that 

a- 2 (y)R(g) - 2A = d(d - 1) (d y In a(y) f , (22) 

where R(g) is the scalar curvature of g a p(x). Therefore, only metrics with constant R(g) 
independent of x can satisfy this equation. They include, in particular, three maximally 
symmetric cases - deSitter (dS d ), plane (IR 6 ' -1 ' 1 ) and Anti-deSitter (AdS ) branes embed- 
ded into the Einstein space with the cosmological constant A. As in the Randall-Sundrum 
case we shall consider the case of anti-deSitter bulk with negative A 

A = —\d{d - l)k 2 , (23) 



so that Eq.(|22|) gives the following profiles of the scale factor a(y) for the foliation of the 



Anti-deSitter spacetime by the deSitter, flat and Anti-deSitter slices 

n 2 



a 2 (y) 



^sinh 2 (%), dS d , 

exp(-2ky), W 1 ' 1 ' 1 , (24) 

^Vcosh 2 (%), AdS d , 
k z 



where the "Hubble" constant TC is given by 

_ \R{g)\ 

n ~ d(d-iy { b) 

3 For reasons of convenience we include the brane tension terms into the bulk part of the action. As far 
as it concerns the Gibbons-Hawking term with the trace of extrinsic curvature, it should also be related 
to the bulk part because it makes by integration by parts the whole Einstcin-Hilbert action quadratic in 
first-order derivatives. 
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Among these three cases the original Randall-Sundrum model [|T^| corresponds to the 
foliation with flat slices. Two flat and empty branes coincide with the pair of such slices 
and serve as boundaries of a smooth piece of Anti-deSitter bulk at y + < y < y_. With 
the orbifold symmetry the piecewise smooth Anti-deSitter solution of vacuum Einstein 
equations satisfies the Israel junction conditions on flat branes when their tensions are 
opposite in sign and equalQ 



/2(d-l)|A| , , 

a+ = - a _=^ 1 - n ' . (26) 

Other cases of curved (deSitter and Anti-deSitter) branes correspond to different ex- 
pressions for their tensions <7j = a±. For Z 2 -orbifold interpolating between the two curved 
branes these tensions, when they satisfy Israel junction conditions, depend on the location 
of branes yi = y± and read 

<7i = -4(d-l)Vn]na(y) , (27) 
where we introduce the notation 

jfl _ / +d y , y = y+, 



1 -o y , y = y-, 



for the inward pointing derivative normal to ±-branes in the background metric fl2T|). 
In order to match with the inward orientation of the normal vector (or location of the 
positive tension brane to the left of the negative tension one) we have to assume that the 
coordinate y in ( f24[) is negative for dS and AdS cases, y + < y_ < 0. This also guarantees 
that deSitter branes are not separated by particle horizon at which a(y) = || [23[ and, 



thus, bound a causally connected bulk. Conversely, relations (|27j) can be viewed as the 
equations which for given bulk curvature parameter k and two different brane tensions of 
opposite signs <r + , <r_ determine the locations of branes at y + and 

Our goal now is to build, along the lines of Sect. 2, the perturbation expansion for 
braneworld action on this background by perturbing the brane configurations like in (|6]) 
and studying the response to it in the bulk ([^). The geometrical picture of the background 
solution above actually shows that a simple scheme of Sect. 2 should be supplemented by a 
number of important modifications due to local diffeomorphism invariance of the theory. 
This invariance manifests itself in many facts. In particular, not all components of the 
bulk metric Gab( x >2/) should be fixed as boundary conditions on the branes, but rather 
the induced metric coefficients g a/ 3, whose number is less than that of GabQ- Moreover, 
the induced metric coefficients are nontrivially related to GAB(x,y) depending on the 

4 For simplicity we dropped in Eq.(|2^) the gravitational coupling constant so that the action has the 
dimensionality of length to the power d — 1 and the brane tensions have the dimensionality of inverse 
length. 

5 The distinguished status of the Gj^-components of the metric that should not be fixed at the bound- 
ary y = const follows from the fact that they serve as Lagrange multipliers of constraints in the 'Hamil- 
tonian' formalism associated with the y-'time' foliation of spacetime. This in turn follows from the local 
gauge invariance of the Einstein theory. 
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embedding of boundaries in the bulk, this embedding being nontrivial after the inclusion 
of gravitational perturbations and depending on the gauge fixation of the latter. 

These properties of the braneworld formalism were considered in much detail in a 
number of papers [14 , 2~T| , [22|| where the location of branes in the bulk was described in the 
Gaussian normal coordinates, G yy = 1, G yfJ- = 0, by an additional scalar field - the radion. 
Here we prefer to proceed in a more gauge invariant manner without explicitly fixing the 
coordinate system in the bulk. The only restriction on the choice of (d + l)-dimensional 
coordinates will be the requirement that the branes are located at constant values of 
the ?/-coordinate^. With this gauge fixation freedom all metric elements are subject to 
nonvanishing perturbations Hab{x, y) which can be decomposed into d-dimensional scalar, 
vector and tensor parts Hab(x,v) = (H yy (x, y), H ya (x, y), H a p(x, y)). The latter two in 
their turn can be decomposed in irreducible scalar, transverse and transverse-traceless 
components with respect to the fixed metric g a /3{x) []• Thus we have 

H yv (x,y) = 2 X (x,y), (29) 
H ay (x,y) = a 2 (y)(v a (x,y) + \y a b(x,y)), \/ p v p {x : y) = 0, (30) 
H a(3 (x,y) = a 2 (y)(>y a p(x,y) + 2 \/( a f/3)(x,y) + 2i/j(x,y)g af3 (x) + 2 \y a \/ e(x,y)), 
\7 a la p{x,y) = g a/3 (x) laf3 {x,y) = 0, XJ a f a {x,y) = 0. (31) 
Perturbations H a p(x 1 y) are generated by the induced metric perturbations h % = h a Jx) 

h %( x ) = a2 (v±) (lap{x) + 24){x)g aP {x) + 2 V(a fp){x) + 2 Va V/3 e (z)) ) ( 32 ) 

as boundary conditions for the linearized Einstein equations in the bulk, ip(x, y±) = ^ ± (x), 
e(x,y±) = e ± (x), f a (x,y±) = f±(x), j a p(x,y±) = 1%( X )- 

We know that the action (|20|) is invariant with respect to (d + l)-dimensional diffeo- 
morphisms. Under their action with the vector field i A {x,y) = (£ a (x, y), £ y (x, y)), whose 
(i-dimensional part has a transverse- longitudinal decomposition in the metric g a /3(x) 

Ca(x,y) = g a p(x)^(x,y) = X a (x,y) + \/ a /i(x,y), \/ a X a = 0, (33) 

irreducible components of (|29|)-(|3T|) transform as 

5X = d y e, (34) 

5b = cT 2 ^ + d y [i, 5v a = d y X a , (35) 

8ip = ip + (d y kia)£ y , 5e = fj,, Sf a (x,y) = X a , (36) 

hoLp{x,y) = 0. (37) 

It is easy to check that the following combinations are invariant under the action of these 
transformations [0] 

. „ , , a~ 2 R(g) „ „ , , a~ 2 

A = Oylp + Ip—r-. r-r—. X Vy l n a j & = Vy e ~ + 1p 



d(d—l)d y \na d y \na' 
V a = d y f a - v a . (3t 



6 At least for small metric perturbations this is always possible, because the foliation by surfaces of 
constant y, Eq.(p^), is regular in patches of the Anti-deSitter spacetime of interest. 

7 Note, that throughout the paper e?-dimensional (Greek) indices are lowered and raised with the help 
of g a p{x) and V« denotes the covariant derivative in the metric g a p(x). 
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The conformal part of the metric perturbation ip(x,y) is not invariant under general 
diffeomorphisms, but for the class of diffeomorphisms leaving the boundaries at constant 
values of the y-coordinate, £ y (x, y±) = 0, boundary values ^(x) are also not transformed. 

One can check that the quadratic part of the action expresses in terms of these invari- 
ants, as it should be for the diffeomorphism invariant functional. Direct calculation shows 
that it decomposes in the sum of decoupled contributions of scalar, vector and graviton 
(transverse-traceless) sectors 



Sg } +1) [H] = S scal [A,B,i[>} + S vect [ V } + S grav [ 7 ] , (39) 



which reac 



UAB^]=2(<i-l) f dV (dA 2 + 2AUB- & HBJ 

B 

^-^J^^-j(n + ^y, (40) 
s Yect [v] = j dvv a (n + ^^j v a , (4i) 

B 

Sgrav[7] = ^ J dV (- d yla p d yl aP + a- 2 laP nj^ - 2 _ a- 2 R{g) laPl ^) . (42) 



Here 



U = g^{x)x 7a x 7p (43) 



is the d- dimensional covariant d'Alembertian. Note that it is defined with respect to the 
y-independent metric g a f3(%) rather than the background metric a 2 {y)g a p{x). Also dV is 
the covariant Riemannian measure in the bulk and dS is the covariant measure on branes 
the latter again defined with respect to the auxiliary metric g a p{x) 



dV = d xdya d (y)y/g(x), (44) 
dS = d d x \J g(x). (45) 



3.1 Scalar and vector sectors 

Interestingly, the calculation of scalar and vector contributions in (|39| ) does not require 
explicit solution of the bulk equations of motion. The reason is that the invariants A, 
B and V Q enter the Lagrangian without y-derivatives and, therefore, the same property 
holds for the variables x, b and v a algebraically entering ([38]). This in turn follows from 
the fact that these variables play the role of Lagrange multipliers of the nondynamical (in 

8 Up to total derivatives in x a directions. Derivatives in y-direction are kept, since they induce terms 
on branes. 
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the extra "time" y) T/A-components of Einstein equations. So, varying the action with 
respect to these Lagrange multipliers we have 



£ S (d+i) 5S (d+1) 

— ^— = -{dyhxa)—^- = -Aa d y^(d-l)(d y \na) (DB + d-A) = 0, 
ox oJ\ 

Jd+l) xa (d+l) 



«sr_ *»sr_ Bto) 



5b 5B v ^ v 7 V d(d 

5S\ d + 1] 5Sg+ 1} rf ( R(g) \ 

In virtue of these equations the vector fl4"I|) and the fritiA; part of the scalar (fjQ|) contribu- 
tions to the action vanishP] and the result reduces to the brane term of fl40|) . 

Note that the scalar sector is completely local and of second order in derivatives. 
However, generally its kinetic terms are not positive definite. With the factors V n In a 
expressed in terms of tensions ( p7|) it takes the form 

S scal [i; ± } = -8(d-l) 2 J2 J dS ^^( n+ JZj)^ ( 47 ) 

which clearly emphasizes its indefiniteness for two tensions of opposite signs, a + > 0, 
o~_ < 0. On the positive tension brane the conformal mode enters with the negative 
(ghost like) sign - exactly as in the Einstein theoryQ, while on the negative tension brane 
its sign is opposite. Only for Anti-deSitter branes with y + < < y- both brane tensions 
can be positive - the case which we will not consider in detail. 

For flat 4-dimensional branes, R(g) — 0, d = 4, this expression reduces to the result 
obtained in f[4|B 



3.2 Graviton sector 

It is important that the tensor (transverse-traceless) sector decouples from the rest of 
variables not only in the action, but already at the level of solving the boundary value 
problem - 7 a a(x) (and only 7 a «(^)) serve as boundary conditions for ^ a p(x,y). That is 
why the situation in this sector literally repeats the technique of Sect. 2. The action in 
the graviton sector (42) after integration by parts reads 



Sgrav[7] = ~ jdV 7 ^(F 7 Q/3 ) + ~ J2j dSad ^ V ^ aP > C 4 * 



9 Up to possible surface terms at the infinity ^"-coordinates in the bulk which we assume vanishing 
because of vanishing sources everywhere except branes. 

10 This does not lead to physical instability because the conformal mode is not dynamically independent 
in view of constraints. 



li 



To compare (47) with the equation (139) of [|l4| one should bear in mind that our conformal modes 



■0 ± are rescaled relative to the conformal modes ^ in jlj], ( y s ± = 2aj_V ± 
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where 

F = a- d d y a% + a' 2 {D - -^—^R(g)^ (49) 

is a second order differential operator that governs the dynamics of gravitons in the 
bulk. Comparison of the equation ( |48f) with ( p~T|) shows that the Wronskian operator 
W± coincides with a±V n on respective branes. Therefore, the solution of the Dirichlet 
problem for transverse-traceless perturbations 

f Fj ap (x,y) = ; 

\ (50) 

in terms of the Dirichlet Green's function of the operator (|49| ) 

FG fl (x, y\x',y') = 5(x,y\x',y') ; 

(51) 

G D {x,y\x',y')\ b± = 0. 

where 5^ d+1 ^(x, y\x', y') = 5 d (x—x')8(y^y')/a d (y)^g(x) is the covariant (d +l)-dimensional 
5-functionQ, reads as 




Substituting it in (f48[) similarly to ( |16D results in the quadratic action in the graviton 




where 

[WGflW]^,!') = a d a d \/ ni V nj G D (x, yi \x', yj ) (54) 

and the primed integration measure is obviously dS' = d d x' ^ g(x'). 

In what follows the operator notations for the kernels of integral operations will be 
more convenient for our purposes. This operator form follows from the fact that the 
differential operator ( [i9"D depends on □ as a parameter, therefore the kernel of its inverse 
can be written down as a function of the □ acting on d- dimensional ^-function 5^(x, x') = 
5 d {x-x')/y/gjx) 

[WG D W]y(D) 5^(x,x') = [WG fl W] i3 (x,x'). (55) 

12 Note that the kernel of the Green's function remains a biscalar with respect to both of its arguments. 
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Thus, finally, after combining together tensor and scalar parts and writing explicitly 
the integration measure, the quadratic part of braneworld effective action looks as follows 



<sf[/i] = ~ IdTxy/gix) V(-7 i (x)[WG s W] 8j (n)7 J (x)+^(x)K iJ (D)^(x)), (56) 




where the notation ip l (x) stands for ip(x,yj) and 

K '^ = « d - 1 ^j4 D + W i T))- (57) 

This matrix operator for the case of d — 4 flat branes coincides via appropriate rescaling 
with the operator K$(d) calculated in [Q by a different procedure. 

At this point it is worth discussing why the sector of transverse-traceless modes in 
(j56|) can be regarded as the graviton sector even for curved branes. The recovery of a 
particle interpretation for the effective action with the essentially nonlocal form factor 
(|55|) is a matter of the low-energy approximation - the expansion of this form factor in 
powers of □ — > 0. In this paper we will not consider this expansion in much detail because 
it was systematically developed in [H] for flat branes. Here we will only sketch a similar 
procedure and its peculiarities due to curvature of a homogeneous brane background. 

To begin with, note that the d- dimensional term □ — 2R(g)/d(d — 1) of the operator 
F, (f49[ ), also enters as a whole the quadratic term 7 Q , / g(D — 2R(g)/d(d — l))7 a/3 of the 
7-perturbation expansion for the rf-dimensional Einstein-Hilbert action 



Seh[9+i}= j d d x ^det(g+ 1 )(R(g+ 1 )-2\) 



(5f 



having a particular value of the cosmological constant 



X = ^2T R{9) > (59) 

which supports the vacuum homogeneous background with a constant scalar curvature 
R(g). This is the reason why in the low-energy approximation the usual Einstein gravity 
theory is recovered on a positive tension brane even when it is not flat. The mechanism of 
this phenomenon is the following. The nonlocal operator ( |55|) of the transverse traceless 
sector parametrically depends on □ — 2R(g)/d(d — 1) as a whole and, therefore, can be 
expanded in powers of this combination. One can show that this expansion starts with 
the first order term (for properly defined graviton fields diagonalising the quadratic form 
in (0), which means that this mode is massless [JIH |23j. The masslessness implies that 
on the maximally symmetric spacetime, R a p^{g) = {g ai ,9fiv - g a u9f3fj,)R(g)/d(d - 1), the 
solutions of the linearized equations of motion for transverse-traceless modes 

D " W^T)) -wW = o (so) 

still have a residual gauge freedom 7 aj a(x) — > j a /3(x) + 2 V(a £p)(x) with the transverse 
vector field, \7c( a = 0, satisfying the equation 2 x/ 13 V(a£/3) = (P + R{.9)/d)i a = 0. The 
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transversality together with this equation guarantee that these transformations preserve 
transverse-traceless nature of 7 a/ g(a;) and the equations of motion fl60|) because (□ — 
2R(g)/d(d — 1)) V(q £a) = 0- Thus, these residual transformations reduce the number of 
physical degrees of freedom for 7 Q/ g by d — 1 to that of the <i-dimensional massless graviton 
d(d — 3)/2. This justifies the interpretation of the transverse-traceless modes as gravitons. 

4 Alternative derivation of the braneworld effective 
action 

There exists an alternative derivation of the effective action suggested in [[13], [TJJ]. It is 
based on the recovery of the action from the effective equations for brane <i-dimensional 
field. These effective equations actually represent the Israel junction conditions on branes 
rewritten entirely in terms of the brane metric itself. In order to see this let us write down 
the metric variation of the fundamental (d +l)-dimensional action (|20|) 

6&[G,(p] = - J d d+1 xVG~(R AB (G)-^R(G)G AB + AG AB J SG AB {X) + 

BxZ 2 

+ J2j ddx ^9 (W - <T#] + i(T"" - gTci)} 5 9lxv {x). (61) 

Here g^ v {x) denotes the induced metric on branes (in contrast with the notations of 
the previous section, where g^ v {x) was reserved for the auxiliary homogeneous metric), 
[K^ u — g^ v K\ is the jump of the extrinsic curvature terms across the brane, and T^ u [x) 
is the corresponding <i-dimensional stress tensor of matter fields on branes 



y/9 5 9, 



fJtU 



X 



Note that the brane part of the ( pT|) gives the Israel junction conditions. 

Suppose that we solved Einstein equations of motion in the bulk subject to fixed brane 
metrics and substituted the obtained solution Gab = Gab[q] in the original action in 
order to find the effective one, S eS [g,<p] = S[G[g],(p]. Then the variation of the latter 



will not contain the bulk term of (plD - because the bulk equations are satisfied - and the 



metric variation of S [g,(p] will be 
8S*[g,<p] 



^g([K^-g^K] + ^(T^-g^c 



<Ti)) ■ (63) 



Here, certainly, the extrinsic curvatures containing the derivatives normal to brane sur- 
faces are the functions of brane metrics. Thus, this equation can be functionally integrated 
to recover S eS [g, </?]. The result ( p3| ) is exact, but expressing its right hand side in terms 
of brane metrics can in general be obtained only within perturbation theory. This proce- 
dure was implemented in |TJj in the quadratic approximation in metric perturbations on 



14 



a flat brane background. To simplify the presentation we briefly repeat this procedure in 
our case. 

We start with the linearized version of the bulk Einstein equations and linearized Israel 
matching conditions on the branes. For transverse-traceless perturbations J a p(x,y) they 
reduce to the linear Neumann boundary value problem 

F7a/3(x,y) = 0, 

(64) 

a d V n -f a/ 3(x,y)\ hi = -~t l ap {x). 
On account of the constraints (H3) in the bulk, the vector part of the equations is trivially 



satisfied, while their scalar part is locally restricted to branes and reads |21 



4(rf-l)^-fn+7^)r'(r) =-/-:,„(,•). (60) 



V ni In a V (d-1] 



From now on we go back to notations of the previous section: F is the operator (49) 
acting in the bulk, V n is the derivative normal to the brane (p8|), t^Jx) is the trans- 
verse traceless part of the rescaled matter stress tensor a d ~ 2 T a p(x) on the i-th brane and 
^-scai = fld_2 9 af3 ( x )Taf3(. x ) (where indices of T^^x) were lowered with the induced metric 
on appropriate brane). Greek indices are raised and lowered with the aid of the auxiliary 
fixed metric g a /3{%) of the homogeneous spacetime, which is conformally related to back- 
ground metrics of all branes (differing only by conformal factors of) which agrees with 
convention of previous section. 

The general solution of Neumann boundary value problem (^4j) is0 

l(x,y) =-nY] I dS ' G N(x,y\x',y i )t i (x'), (66) 



2^ 

1 b 

where Gjq{x,y,x\yi) is the Neumann Green's function 

FG N (x, y\x', y') = 5^ d+1 \x, y\x', y'), 

V n G N (x,y\x',y')\ b _ = 0. 



(67) 



Effective dynamics of the brane metric allows one to express the latter solely in terms 
of brane matter sources. In the transverse-traceless sector this reduces to restricting the 
general bulk solution to branes, i.e. by putting in (|66| ) y = yj. In terms of the operator 
notation for the kernel of the Neumann Green's function (similar to 



G N {x, yj \x', yi ) = G^(D)S^{x,x'), (68) 

this reads as 



13 As in the case of the Dirichlet boundary value problem above in what follows we drop for brevity 
Greek indices in linear and quadratic combinations of 7 Q( a(a;,y) and t^p{x). 
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There are infinitely many actions that generate by variational procedure the last equa- 
tion. Fortunately, there is a hint coming from Eq. (|63| ) that allows one to fix this ambiguity 
- the rescaled matter stress tensor enters the variational derivative of the braneworld ac- 
tion with respect to induced metric a 2 g a p(x) with the algebraic coefficient a~ 2 y/g(x). 
Therefore, by rewriting (|69| ) and ( |65D as 

( G n)v ( D )^» + \yW)U(x) = 0, 

j 

^gjx) + ^gjx)tr\x) = 0, (70) 

j 

where Ky(d) is defined by (|57|), and (G^ 1 ) .^.(n) is the inverse of G^(D) 

£ (G^ 1 ) <,(□)(#(□) = 6}, (71) 
j 

one easily finds the action generating ( [70"D 




b 



The source term with i J (x) and its trace here is obviously a linear in metric perturbation 
part of the action of matter field, while the quadratic in 7 part has a purely gravitational 
nature and, thus, should necessarily coincide with the action (|56|) obtained within the 
Dirichlet boundary conditions. This implies the equality of two operators arising in two 
different boundary value problems - the Dirichlet and Neumann one 

[WG D W]^n) = -(G w %(D). (73) 

5 Duality of boundary value problems 

To prove the duality relation ([73]) let us, first, remind the notations of Sect. 2. We 
consider At - the bulk spacetime with branes/boundaries dAif dAi = [J i dAii. Let X 
be the coordinates on the bulk M. and x - the coordinates on the branes dAii, so that 
the embedding of dM.i into M. reads as X = Xi(x). Introduce H - the field in the bulk 
and denote boundary values of the bulk field on dAii by h\ h l (x) = H(Xi(x)). Let F 
be some nondegenerate self-adjoint differential operator of the second order in derivatives 
acting on H(X). Let it satisfy the Wronskian relation (|T0| ) with the first order Wronskian 
operator W. 

Let H(X) satisfy the Neumann problem with the inhomogeneous boundary conditions 
with some given sources jk{x) on the boundary 

FH(X) = ; 

(74) 

WH(X)\ x=Xk{x) = Jk (x). 
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The solution of this problem in terms of the Neumann Green's function of the operator 
F, satisfying^ 

f FG N (X,X') = 6(X,X') 1 

(75) 

{ WG N (X,X>)\ x&aM =0, 



reads 



H(X) = J2 I dS'G N (X,X k (x')) Jk (x'). 



(76) 



l! 8M k 



Its restriction to the boundary yields the relation between the boundary fields h (x) 
and their conjugates - sources j k (x) 



h\x) = Y. I dS'G l *(x,x')j k (x') (77) 



k 9M k 



in terms of the new boundary operator kernel 

G%(x,x') = G N (X l (x),X k (x')) (7* 



On the other hand the same field ( ffq) serves as the solution of the Dirichlet problem 
(H) with the boundary data H(X)\ x =x{x) = h(x) related to the sources j(x) according 



to ([77| ). Therefore, in terms of the Dirichlet Green's function flT2|) H{X) has the form 
(|l3l). Acting on (|T3| ) on the boundaries by the Wronskian and taking into account that 
WH(X)\ gM = ji(x) one finally obtains the alternative to (|77|) relation between h(x) and 



) = ~J2 / dS'[WG D W} kl (x,x')h\x'), (79) 



i 



dM 



where the kernel [WG D W] H (x, x') is defined by Eq.(|19l). 

Comparing (fT"7|) and (|79|) one concludes that G l ^(x,x') and —[WGDW} k i(x,x') are 
the kernels of the inverse operations on dAi 

-J2 J dS'{WG D W} H (x,x")G%(x,x') = 5 l k 5{x,x'). (80) 

This confirms the operator relation (|73|) . Below we consider two examples demonstrating 
the efficiency of this result. 

14 Wc define the delta-function 5{X,X') by the relation / dV H{X) S(X,X') = H(X'). 

M 
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5.1 Simple example 

Consider the simplest model with the operator on the one-dimensional space of the 
coordinate y, y\ < y < y2, 

" ^ I nr. (81) 



dy 2 

and introduce its Wronskian operator acting on boundaries 

+d v , V = Vi, 
-d y , y = y 2 . 

The Dirichlet and Neumann Green's functions of this operator read 



\2) 



smm(y'-y 2 )smm(y-yi) , , , 

G D {y,y) = : — - f v &{y -y) + {y <->y), 83 

msinm(i/2- y-y) 
n ( >\ cosm (y' -2/2) cos m(y-yi) , 

G N (y,y) = : r 6(y -y) + {y ++y), 84 

m sin m(y2 — yi) 

so that the operators (|T9|) and (WW) turn out to be 2 x 2-matrices 



-> «- m f-cosm(y 2 -yi) 1 

-[V„G D V n ] i3 - = ? r , (85) 

smm(y 2 -yi) \ i / > 

V»z yiy y ^ _ cosm ^ 2 _y 1 ^ 

/cosm(y 2 -yi) 1 

= — 7 V I • (86) 

msmm[y 2 —yi) \ r \ 

w a J y 1 cosm(y 2 — yi) 

This is a matter of direct verification that these matrices are inverse to one another, 
because their determinants equal 

det G N ij = (det[WG D W]^) _1 = --^ (87) 

5.2 Example: two-brane Randall-Sundrum model 

In this section we illustrate the relation (173) for the Randall-Sundrum model with two 
flat branes. As in Sect. 3 we work in coordinate system X A = (x a , y) in which background 
metric has the form ( pl~l) with flat spacetime metric g a p(x) = i] a p. Two flat branes with 
tensions cr + = — cr_ are located respectively at y = y + and y = y_. In such a coordinate 
system 

F = a- d d y a% + a~ 2 n ; □ = if* V. V/3 , (88) 
\x)r] a/3 (x) = d , a(y) = exp(-2%) , 
.Avy _ / +d y , y = y+, 



V n = n A V A ={ y " y+ ' (89) 

1 -9 tf , y = y- 
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For the sake of brevity we introduce the new operator-valued variable depending on 
the coordinate y 



x = x(y, U) 



ka(y) ' 

x + = x(y+, □) , x_ = □) , (90) 

and the following two-point combinations of Bessel functions 

D(x, x ) = - (ld /2 (x) J d/2 (x') - Krf/a^) J d/2 (x)) , 

E(x, x') = - (Y d/2 (x)J d /2-i(x') ~ Yd/2-i(x)Jd/2(xj) , 

C(X, X') = I (y d/2 _l(x) J d/2 -l(x?) - Y d/ 2-l{^)J d /2-l{^)) ■ (91) 

Then the explicit expressions for Dirichlet and Neumann Green functions take the form 



3 



n , n n , D(V, x_)D(x, x + ) Ojy'-y) + D(V, x + )D(x, x_) %-j/Q 

GD(2/ ' y |D) = a(2/) ,/ 2a(?/r / 2fcD(x _ )X+) > ( 92 ) 

n , E(x',x_)E(x,x + ) %'-?/) + E(x / ,x + )E(x,x_)g(y-y / ) 



a(|/) d / 2 a(|/O d / 2 A;C( 



x_, x 



After acting on the Dirichlet Green's function with normal derivatives fl89|) and re- 
stricting them to branes (y,y f = y +) yJ) one obtains the 2 x 2-matrix valued operators - 
the kernel of the braneworld effective action in the Dirichlet setup 

[WG,W],(D) = k I | , (94) 

D;X -' X+) \ afa d J 2 alx_E(x + ,x_ 



'a^[x + E(x_, x + ) a+a_ 



and the inverse of the kernel in the Neumann setting (|7 

/, s-i /— a7 d x_E(x + , x_) a, d ^ 2 a_ d ^ 2 

G « (n) = ^) |. (95) 

\ a + a_ — a_ x + E(x_,x + ; 



Here we used the properties of Bessel functions p5[ 

^(x d/2 D(x' ; x)) = x d/2 E(x', x) , E(x, x) = -x" 1 . (96) 

The calculation of the 2 x 2-determinants of these matrices can be easily done by using 
a simple relation 

E(x', x')E(x, x) - E(x', x)E(x, x') = D(x, x')C(x, x') (97) 
from which it follows that 

detGg(D) = (detfWG^^n))" 1 = - ^^f' • ( 98 ) 

Then it is a matter of direct check that the matrices (|94]) and (|95l) are inverse to one 
another. 
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6 Conclusions 



Thus we obtained the braneworld effective action in the weak field (quadratic) ap- 
proximation for models of curved branes with the deSitter and Anti-deSitter geome- 
tries, embedded in the Anti-deSitter bulk. The calculation was done within two dif- 
ferent schemes resorting to Dirichlet and Neumann boundary value problems, and their 
equivalence was shown on the basis of a special relation between the boundary operators 
associated respectively with the Dirichlet and Neumann Green's functions of the theory. 
Although this relation was used here only in a braneworld context, it seems to have a much 
wider scope of implications in various models relating volume (bulk) and surface phenom- 
ena, like graviton localization [O, AdS/CFT correspondence and holography principle 



27, 28, 29 



Regarding the braneworld action algorithms for curved branes - this result can be 



important for the brane models of inflation scenario [11, 13, 14 1. One of the motivations 



for considering deSitter type branes is that they model the inflation stage in the dynamics 
of the brane Universe embedded in the multidimensional Anti-deSitter bulk. The idea of 
the radion field playing the role of the inflaton was suggested in |0|, where inflation was 



induced by the effectively repulsive force between the branes diverging in the course of 
inflationary stage. This repulsive force arises in virtue of the detuning of brane tensions 
from their Randall- Sundrum values corresponding to flat branes. However, the anal- 



ysis of this detuning in |L3], [L4] was done by directly extrapolating the flat branes results 
to curved branes. Now, with the knowledge of the answer ([56]) the effects of curvature can 
be studied at a rigorous quantitative level. This and, in particular, the effects of brane 
curvature on the low-energy limit of the nonlocal form factor (54)- (55) will be considered 
elsewhere. 

Another possible generalization of the obtained results concerns going beyond the 
limitation of conformally equivalent branes. Note that everywhere in the constructions 
above both branes are "parallel" surfaces simultaneously of the flat, deSitter or Anti- 
deSitter geometry. It would be interesting to generalize calculations to coexisting branes 
of different types (different signs of curvature), because this might lead to the dynamical 
description of intersecting curved branes very popular in recent years [0, EH], Q . 

Finally, the algorithm (|56|) suggests a natural generalization to multibrane cases when the 
range of index i goes beyond two possible values. This generalization particularly implies 
that a brane can serve as a boundary of more than two smooth bulks glued to it. Such 
a construction of branching bulks can be shown to have explicit Randall-Sundrum type 
solutions of Einstein equations with appropriately generalized Israel junction conditions 



| 35fl . It forms a kind of spacetime bulk-brane foam resembling selfreproducing inflationary 



Universe [B6| and suggesting new interesting facets in the cosmological constant problem. 
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